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Abstract. 

We show theoretically that two evanesccntly coupled x^^^ second harmonic 
generators inside a Fabry-Perot cavity provide a tunable source of quadrature squeezed 
light, Einstein-Podolsky- Rosen correlations and quantum entanglement. Unlike 
systems using coupled downconvcrters, second harmonic generation has no oscillation 
threshhold, so that the entangled fields become macroscopically occupied as soon as 
the pumping fields are turned on. This system also gives two frequencies at which 
the entangled fields can have macroscopic intensity. We show how the entanglement 
properties can be controlled by adjusting the pumping, the coupling strengths and the 
cavity detunings. 
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1. Introduction 

The system of two evanescently coupled x^'^^ processes operating in the second harmonic 
generation (SHG) regime inside a pumped Fabry-Perot cavity has been theoretically 
analysed by Bache et al [1], in terms of suppression of number difference fluctuations 
between the output modes. In this work the name quantum optical dimer was given 
to this device. Although the calculations performed by Bache et al were sufficient to 
show that strong intensity correlations should exist between the outputs, they did not 
calculate any of the phase-dependent correlations which are necessary to demonstrate 
the existence of continuous variable bipartite quantum entanglement. In later works, 
Olsen and Drummond [2] and Olivier and Olsen [3] have shown theoretically that a 
related device, operating in the downconversion regime of optical parametric oscillation, 
can provide a robust source of continuous variable entanglement and be used to 
demonstrate the Einstein-Podolsky-Rosen (EPR) paradox in both the above and below 
threshold regimes. In this work we will examine the quantum optical dimer proposal of 
Bache et al in terms of the production of entanglement and EPR states. In principle this 
sytem may have all the operational advantages of the coupled downconverter proposals, 
along with the property that all output modes are always macroscopically excited in 
second harmonic generation since there is no oscillation threshold. 

Generally, the device we are proposing may be considered as either a single 
nonlinear crystal pumped by two spatially separated lasers, or two waveguides with 
a x^^^ component. We calculate phase-dependent correlations between the outputs of 
the cavity, evaluating entanglement criteria due to Duan et al Simon [5], the 

lograithmic negativity j6], as well as the Reid-Drummond EPR correlations [3, [H [9] 
which can also be used to prove entanglement. Generically, our system is related to 
the nonlinear coupler, which is the name given to a system of two coupled waveguides 
without an optical cavity by Perina et al [10]. The device generally consists of two 
parallel optical waveguides which are coupled by an evanescent overlap of the guided 
modes. The quantum statistical properties of this device when the nonlinearity is of the 

type have been theoretically investigated, predicting energy transfer between the 
waveguides [11] and the generation of correlated squeezing [12] . When operated inside an 
optical cavity, entanglement between the output modes has been predicted [I3] . Here we 
show that the system with x^'^^ nonlinearity, and operating in the upconversion regime, is 
also potentially an easily tunable source of single-mode squeezing and entangled states 
for both the low frequency fundamental modes and the high frequency up-converted 
harmonic modes. The spatial separation of the output modes means that they do not 
have to be separated by optical devices before measurements can be made, along with 
the unavoidable losses which would result from this procedure. We will show that the 
correlations are tunable by controlling some of the operational degrees of freedom of 
the system, including the evanescent couplings between the two waveguides, the cavity 
damping rates, the input powers and the cavity detunings. 
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2. Entanglement measures 

Before we begin to analyse the actual physical system, we will define and outline 
the measures which we will use below to demonstrate continuous-variable bipartite 
entanglement. Entanglement is a property of quantum mechanics which is related to 
the inseparability of the combined density matrix of a system into density matrices for 
its subsystems. In the present situation, we will be interested interested in continuous 
variable bipartite entanglement between the output modes from each side of the dimer, 
which we shall label as 1 and 2. It is firstly necessary to give the definition of the optical 
quadratures we will use, as the exact form of any inequalities depends on this. For the 
description of this system, we require four intracavity bosonic annihilation operators, 
Oi, 0,2, bi, b2, where the cij are for the low frequency (fundamental) modes and the bj 
are for the high frequency (harmonic) modes. We define the quadrature operators at 
the phase angle 9 as 

Kj = + (1) 

so that j, -^ff^^^] ~ 2z5jfc, where 6jk is the Kronecker delta, and similarly for the 
high frequency quadratures. The Heisenberg uncertainty principle (HUP) then requires 
V {X^ j)V {X^'^'^^'^) > 16jk- In the interests of notational simplicity, we will label the 

quadrature -^f ^'^^^ as Y^i^- 

The EPR paradox stems from a famous paper published in 1935 [H], which 
used a gedanken experiment with particles which were entangled in position and 
momentum to show that local realism was not consistent with the completeness 
of quantum mechanical theory. A direct and feasible demonstration of the EPR 
paradox with continuous variables was first suggested using nondegenerate parametric 
amplification [151 CSl [II]- This was possible because the optical quadrature phase 
amplitudes used in these proposals have the same mathematical properties as the 
position and momentum originally used by EPR. Even though the correlations between 
these are not perfect, they are still entangled sufficiently to allow for an inferred violation 
of the uncertainty principle, which is equivalent to the EPR paradox [3 O [9]. An 
experimental demonstration of this proposal by Ou et al soon followed, showing a clear 
agreement with quantum theory [18]. A recent theoretical proposal has examined a 
demonstration with the atomic field quadratures of massive particles, using the coherent 
dissociation of molecular Bose-Einstein condensates [T9] . 

Although the concept of entanglement is required to formulate the EPR paradox, 
for completeness we will outline the details of a proof given by Reid [20] that seeming 
violations of a Heisenberg uncertainty principal (HUP) as in the experiment of Ou et al 
[T8] are automatically demonstrations of entanglement. We begin by assuming that a 
given system is bipartite separable and divide it into two subsystems A and B. We now 
consider observables x"^ and of subsystem A, obeying the Heisenberg uncertainty 
principal with V{x^)V{ii^) > 1. We now introduce Vinj^XA) as the measured error 
in the prediction for the outcome of a measurement x^ at A, based on a result at 
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B, and similarly for VinfivA)- The first task is to show that separability always 
demands Vinf{x'^)Vinf{y'^) > 1, so that violation of this inequality requires inseparability 
and hence entanglement of A and B. The conditional probability of result for a 
measurement of at A given a simultaneous measurement of x^ at B with result xf 
is 

nA^f) = ^^f^, (2) 

where, assuming separability 

P(x^, xf) = P{r)P{xf\r)P{xy), (3) 

r 

with a separable density matrix being written as 

p = J2P(r)pfp?. (4) 

r 

In this case, \x^) and \x^) are eigenstates of x^ and x'^, with P(x'^|r) = {x^\p^\x^) 
and P(x'^|r) = {x^\p^\x^) . The mean of this conditional distribution is 



X 



_ ZrPir)P{xf\r){x^)r 
P{xf) 

where {x^)r = ^^.a P{x^\r). The variance, Vi{x), of the distribution P(a:'^|xf ) is then 

P(r)P(xf |r)^^P(x^|r)(x^ - /i^ 

For each r, the mean square deviation, 'Y^^a P{x^\r){x^ — dY , is minimised by the choice 
d = {x^)r, so that for the choice d = pi, 



where is the variance of P{x^\r). We may also define a measured error, 

Vinf,est{x^), resulting from linear inference, which will not be better than that based 
on knowledge of the conditional probabilities, so that 



X 



P{X, 



^P(r)Ki 



xf 



X 
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We also note that Vinf{y^) > J2r^(''^)^r{y^), where P{y^\r) = {y^\Pr\y^) and \y^) is 
the eigenstate of y^. The Cauchy-Schwartz inequahty then imphes that 

V.nf{x^)V,^f{y^) > Y,nr)Vr{x^) x Y.P{r)Vr{y^) 

(9) 



r 



For any physical p,., the HUP requires that Vr{x^)Vr{y^) > 1. For a separable state it 
is therefore required that 

Vinf{x^)V,nf{y^) > 1. (10) 

As the measured errors are always at least as large as the inferred errors, this means 
that an observation of Vinf^est{x^)Vinf,est{y^) < 1 is sufficient to prove inseparability 
and hence bipartite entanglement. We note that this violation is a sufficient but not a 
necessary condition, so that continuous variable bipartite entanglement may be present 
which is not detected by this means. An example of such a situation, which can occur 
for mixed states, is given by Bowen et al [21]. We also note here that the EPR paradox 
signifies a stronger form of entanglement, as has been discussed recently by Wiseman et 
al in the context of steering |22j . 

In practice any EPR measurement is usually done by defining orthogonal 
quadratures via a minimisation of the errors in a linear inference process [23]. The 
apparent violation of the HUP for these quadratures then signals entanglement between 
the modes of the system. We note here that it has also been shown by Tan [21] that 
the existence of two orthogonal quadratures, the product of whose variances violates 
the limits set by the HUP, provides evidence of entanglement. Tan demonstrated 
this in the context of teleportation, with the outputs from a nondegenerate OPA 
mixed on a beamsplitter. This procedure has been extended to the case of tripartite 
entanglement [25] and its extension to larger numbers of modes is straightforward, even 
though the number of inequalities to be violated increases due to the different classes of 
inseparability which then exist [26] . 

The second of the entanglement measures is due to Duan et al [4] and also Simon [5], 
who developed inseparability criteria which are necessary and sufficient for Gaussian 
states, and sufficient in general. These criteria have recently been shown to be special 
cases of an infinite series of inequalities based on the non-negativity of determinants 
of matrices constructed from certain combinations of operator moments [27]. In the 
general case, we may define combined quadrature operators similarly to Duan as 

xi = \b\xf ± 

Yl = ± If/, (11) 

where b is an arbitrary non-zero real number. It may be shown that, for separable states, 
V{Xi) + V{Y^) >2(b' + ^Y (12) 
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with any violation of this inequahty therefore demonstrating the presence of bipartite 
entanglement. In what follows, we will choose fe = 1 so that the lower bound of the 
inequality is 4. While this is not the optimal choice for the general case, it is the 
appropriate choice for this system, due to the symmetry between modes 1 and 2. 

The third measure which we will apply is the logarithmic negativity, proposed by 
Vidal and Werner as a computable measure of entanglement, as opposed to others which 
can be difficult to calculate [6]. We note here that this measure is defined for Gaussian 
states, to which we are also limited here due to the linearisation process which we will 
be using. We first define the system covariance matrix as 



where 



and 



Defining 



C 



Ci 

C21 



C12 
C2 



V{Yf,Xf: 



V{X. 
ViYj 



Yf 



V{X, 
V{Y, 



■ Y^ 



^ = ^ (det Ci - det C12) - V(det C2 - det €12^- det C 
the logarithmic negativity is then defined as 



(13) 



(14) 



(15) 



(16) 



HO 



-loga^ if e < 1 
otherwise. 



Any non-zero value of ^(0 is then an indication that the two modes are entangled. An 
interesting feature of this measure is that it has no dependence on quadrature angle, 
but becomes a function of frequency only. This shows that the logarithmic negativity is 
useful for demonstrating that continuous variable entanglement exists in a given system, 
but does not tell us at which quadrature angles the system may exhibit the necessary 
properties for uses such as teleportation. 



3. The system and equations of motion 

The physical device we wish to examine is the same as that described in reference [1]. 
As this device has been described there in detail, we will give a briefer description 
of the essential features here. The system consists of two coupled nonlinear x^'^^ 
waveguides inside a driven optical cavity, which may utilise integrated Bragg refiection 
for compactness. Each waveguide supports two resonant or near resonant modes at 
frequencies Ua (fundamental) and Ub (harmonic), where 2uJa — ujf,. The lower frequency 
modes at Ua are driven coherently with external laser fields, while the nonlinear 
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interaction within the waveguides produces second harmonic photons with frequency 
Uh. We assume that only the cavity modes at these two frequencies are important 
and that there is perfect phase matching inside the media. The two waveguides 
are evanescently coupled. We will be interested in the phase-dependent correlations 
necessary for an unambiguous demonstration of entanglement and the EPR paradox, 
rather than intensity correlations considered in reference [1]. 
The effective Hamiltonian for the system can be written as 

where the nonlinear interactions with the x*^^"* media are described by 



In the above k denotes the effective nonlinearity of the waveguides (we assume that 
the two are equal), and dk, bk are the bosonic annihilation operators for quanta at the 
frequencies Ua, uJb within the nonlinear medium k (= 1, 2). The coupling by evanescent 
waves is described by 

bibl + b\k], (19) 



'^couple 



did\ + a[d2 



hJh 



where the Jk are the coupling parameters at the two frequencies, as described in reference 
[1], where it is stated that the lower frequency coupling, J^, is generally stronger than 
the higher frequency coupling, J^, and also that values of Ja as high as 50 times the 
lower frequency cavity loss rate may be physically reasonable. The cavity pumping is 
described by 



T~(pump 



eia[ — e*ai + e2dl — 62^2 



(20) 



where the ek represent pump fields which we will describe classically. Finally, the cavity 
damping is described by 



7ir 



-^T.{^'al + rfbl)+h.c., (21) 

k=l 

where the represent bath operators at the two frequencies and we have made the 
usual zero temperature and Markov approximations for the reservoirs. 

With the standard methods [28], and using the operator/c- number correspondences 
{dj aj, bj the Hamiltonian can be mapped onto a Fokker-Planck equation for 

the Glauber- Sudarshan P-distribution [291 EH]- However, as the diffusion matrix of 
this Fokker-Planck equation is not positive-definite, it cannot be mapped onto a set of 
stochastic differential equations. Hence we will use the positive-P representation [31] 
which, by doubling the dimensionality of the phase-space, allows a Fokker-Planck 
equation with a positive-definite diffusion matrix to be found and thus a mapping 
onto stochastic differential equations. Making the correspondence between the set of 
operators (a^, aj, foj) (j = 1,2) and the set of c-number variables {aj^a'j' , Pj, we 
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find the following set of equations, 



^ = ei - (7a + iAa)ai + nafPi + iJaa2 + \f^\ r]i{t), 



+ 



dt 
da2 



dat 



£2 - (7a + iAa)a2 + ^2 + iJa^l + Vsit) , 

'- ^2 - (7a - iAa)a^ + Ka2(3t - '^JaOit + \J uPt Vi{t)^ 



dt 

_4 
dt 

= -{lb + iAfe)/yi - + iJbP2, 

dpi 



(22) 



— - -(76 - iA,)/3+ - 2 «i - i 

-T- = -(7b + iAfe)/y2 - -«2 + ^-JbPl, 
dt 2 

^ = -(7. - iA,)/5+ - ^at ' - U,Pt, 

where the 7^ represent the cavity damping rates at each frequency. We have also added 
cavity detunings Aa^b from the two resonances, so that for a pump laser at frequency 
ul, we have Aa = uja — ojl and A;, = c^;;, — 2ujl- Below, in section [6], we will investigate 
detuning effects in greater detail. The real Gaussian noise terms have the correlations 
rij{t) = and 'Uj{t)rik{t') = 6jkS(t — t'). Note that, due to the independence of the noise 
sources, {Pk) and {P'^) are not complex conjugate pairs, except in the mean over 
a large number of stochastic integrations of the above equations. These equations allow 
us to calculate the expectation values of any desired time-normally ordered operator 
moments as classical averages, exactly as required to calculate spectral correlations. 



4. Linearised analysis 



In an operating regions where it is valid, a linearised fluctuation analysis provides a 
simple way of calculating both intracavity and output spectra of the system [32], [33], 
by treating it as an Ornstein-Uhlenbeck process [31]. To perform this analysis we first 
divide the variables of fl22l) into a steady-state mean value and a fluctuation part, e.g. 
«! — > af^ + 6ai and so on for the other variables. We find the steady state solutions 
by solving the equations (l22l) without the noise terms (note that in this section we will 
treat all fields as being at resonance), and write the equations for the fluctuation vector 
Sx = [Sai, 6af, Sa2, Sa2, Sj3i, SP^, SP2, ^PtV^ ^^^^ order in these fluctuations, as 



d6x = -A6x dt + BdW, 
where the drift matrix is 



A 



Aaa 
Aba 



-A 
Ab 



ba 



(23) 



(24) 
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Figure 1. The single-mode output variances for Xa and Xb, for ei = £2 = 0.8ec, 
Aa ~ Ab ~ 0, -fa = 7b, and different values of Ja and Jf,. The solid lines are for the 
uncoupled case, with Ja = Jb ~ 0. The dash-dotted lines are for Ja = Jb ~ Ja, with 
the optimal quadrature angles being 9 = 18° (fundamental) and 13° (harmonic). The 
dashed line is for Ja = 2ja ~ "^Jb and 9 = 46° (fundamental) and 74° (harmonic). A 
value of less than one indicates squeezing. All values plotted here and in subsequent 
graphics are dimensionless. 



with 



la + i^a 




-iJa 







la - i^a 





iJa 


-iJa 





la + i'^a 







iJa 




la - i^a 



(25) 



and 



A, 



ba 

















Kaf * 




























Kttf * 



(26) 



and 



A, 



bb 



lb 





-iJb 








lb 



lb 



-iJb 






7b 





iJb 


-iAh 



(27) 



In 





iJb 

dW is a vector of real Wiener increments, and the matrix B 
is zero except for the first four diagonal elements, which are respectively 
■\/k^, -^/k^^P*, \J k/?!"*, \J The essential conditions for this expansion to be 
valid are that moments of the fluctuations be smaller than the equivalent moments of 
the mean values, and that the fluctuations stay small. In the case of uncoupled SHG, it 
is well known that there is a critical operating point above which this condition does not 
hold and the system enters a self-pulsing regime [35], [361 Elj- This point is easily found 
by examination of the eigenvalues of the equivalent fluctuation drift matrix for that 
system, and this procedure is also valid in the present case. The fluctuations will not 
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tend to grow as long as none of the eigenvalues of the matrix A develop a negative real 
part. At the point at which this happens the linearised fluctuation analysis is no longer 
valid, as the fluctuations can then grow exponentially and the necessary conditions 
for linearisation are no longer fulfilled. In this work we will only be interested in a 
region where linearisation is valid and will restrict our analyses to below the self-pulsing 
threshold, found for 

e = e, = ^2:^V276(7a + 7.) (28) 
in the case of a single uncoupled cavity. 




Figure 2. The spectral correlations 5™*(X_) + S'^"*(y+) of the fundamental (left) and 
harmonic (right), for ei = £2 = 0.8ec, = At, = and different values of Ja and Ji,. 
The solid lines arc for Ja = Jb = Ja^ with the optimal quadrature angles being 9 = 63° 
(fundamental) and 58° (harmonic). The dash-dotted lines are for = 27^ = 2Jf,, at 
6 = 91° (fundamental) and 114° (harmonic), while the dashed lines are for Ja = 27^ 
and Jf, = 7a/2 at = 179° (fundamental) and 119° (harmonic). A value of less than 
4 indicates bipartite entanglement. 



To examine the stability of the system, we first need to find the steady state 
solutions for the amplitudes, by solving for the steady state of ( |22|) with the noise 
terms dropped. For this system in the general case, we find that it is convenient to 
solve for these steady states numerically, using a Runge-Kutta algorithm to integrate 
the system of deterministic equations until the solutions are well into the steady-state 
regime. The eigenvalues of the drift matrix are then also found numerically. Using 
the steady-state solutions, we may then calculate any desired time normally-ordered 
spectral correlations inside the cavity using the simple formula 

S{uj) = {A + iujiy^ BB^ (A^ - icul) , (29) 

after which we use the standard input-output relations [33] to relate these to quantities 
which may be measured outside the cavity. For example, the spectral output variance 
of the quadrature Xf ± X2 will be denoted Sg^^{X±), with whether it refers to the 
fundamental or harmonic being made obvious by the context. In what follows we will 
use values oi k = 0.01 and 7a = 1 while varying other parameters. 
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Figure 3. The spectral correlations S°^j{X)S°^j{Y) of both frequencies, for ei = 
£2 = 0.8ec, Aq = Ab = and different values of Ja and Jf,. The solid lines are for 
Ja = Jb = la, with the optimal quadrature angles being 6 ~ 174° (fundamental) and 
45° (harmonic). The dash-dotted lines are for Ja = = 2 Jb a.t 9 = 46° (fundamental) 
and 24° (harmonic), while the dashed lines are for Ja = 27a and Jb = 7a/2 at 6* = 44° 
(fundamental) and 29° (harmonic). A value of less than 1 illustrates the EPR paradox 
and hence bipartite entanglement. We see that this measure is markedly less sensitive 
to entanglement at the lower frequency than the correlations of figure [21 

5. Doubly resonant cavity 

In this section we will give results for the case where all four intracavity fields are at 
resonance. The first correlations we show, in figure [1] are the minimum quadrature 
variances for each mode for different values of the coupling strengths. We see that the 
quadrature angle of minimum noise is changed by comparison with the resonant case, 
and that the harmonic exhibits no single-mode squeezing for the coupling strengths 
shown. The single-mode squeezing in the fundamental is also noticeably degraded 
from the uncoupled configuration, and both show excess noise at some frequencies. 
This is common with entangled systems, where the modes considered individually will 
show excess noise and the entanglement manifests itself in the correlations of this noise 
between the two modes. The fact that the minimum values of the correlations are no 
longer found at either ^ = 0or7r/2isa consequence of the evanescent coupling, and 
has previously been seen in the coupled downconversion configuration [21 [3] . 

In the next three figures we examine the Duan and Simon criteria, the EPR criteria 
and the logarithmic negativity as the couplings are changed. In figure [2] we show 
the minima of the correlations at the quadrature angles for which these are found. 
We see that the angles are quite different for the fundamental and harmonic for the 
same parameters and that the degree of violation of the inequalities varies markedly 
with coupling strengths. The fact that the quadrature angles change should not be a 
problem experimentally, as homodyne detection techniques generally scan through all 
angles. When we examine figure [3], we see that the EPR criteria do not provide a very 
sensitive measure of the presence of entanglement for the fundamental modes but are 
better for the harmonic, although the entanglement between the two harmonic modes 
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0.8 
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Figure 4. The logarithmic negativity, for ei — £2 — O-Scc, = = and 

different values of Jq and Ji,. The left plot shows the values for the low frequency 
modes and the right plot shows them for the high frequency modes. The solid lines 
are for Ja ^ Jb = Ja, the dash-dotted lines are for Ja = 2ja ~ 2 J;,, and the dashed 
lines are for J a — "i^a and Jf, = 7a /2. 



for Ja = Jb = "Ja is completely missed. This is not in contradiction with the results 
shown in figure [2] because the EPR criteria are merely sufficient but not necessary to 
demonstrate inseparability of the modes. It is also apparent that they give their best 
results at different quadrature angles to the Duan and Simon criteria. 

The logarithmic negativity is perhaps the most sensitive indicator of inseparability 
for this system, although it does not tell us the quadrature angle at which maximum 
violations of the other inequalities can be found. In ffgure HI we see, for example, that 
it is much more sensitive at finding entanglement between the fundamental modes than 
the other two methods in the case where Ja = Jb = la- We note here that the Duan 
and Simon criteria showed very little violation of the inequality in this case, but we 
have shown only the quadrature angle of maximum violation. The frequency range over 
which the logarithmic negativity is positive tells us that violations can be found for other 
quadrature angles as the frequency changes, although given the fact that the maximum 
violation of the Duan and Simon inequality is very small here, it is unlikely that this 
particular combination of parameters would be a good operational choice. 

6. Detuning the cavity 

Often in optical systems the best performance is found when the cavity is resonant 
for the different modes involved in the interactions, with any detuning worsening the 
quantum features such as squeezing and changing the quadrature angles at which these 
are seen [38]. In the present case we find that detuning the cavity by the appropriate 
amount from the two frequencies allows for some simplification of the theoretical analysis 
and can actually improve some quantum correlations. In the special case where = Ja 
and Ab = Jb the best quadrature angles remain fixed as the pumping is varied. If we also 
set ei = e2 = e, the equations become somewhat simplified and we are able to progress 
further analytically. In this case we define the new steady-state variables, a± = aii: a2 
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and (3± = (3i±(32 and solve the classical equations for these. Due to the symmetry of the 
system, we may assume that a_ = /5_ = 0, and that the variables are real. Numerical 
analysis bears out these assumptions. We then find that 



— - — a 



47b 



with being the real solution of the cubic equation 



Setting 



we find 



K 

8^1 



X 



■a+ + laa+ - 2e = 0. 



2X ^lalh 



1/3 



(30) 



(31) 



(32) 



(33) 



3V3;^- 

We now need to solve for the fluctuations in the new variables. Setting ai + 
02 = a+ + 5«+ etc, we write the equations of motion for the variables 6x± = 
[Sa+, 6a~\^, 6a-, 6a1l, Spj^, 6P-, SP^]'^ in matrix form as 



where 



A 



pm 



6x± 




Apm6x±dt + B±dW, 








la 


-f/3- 


f 





—-0* 
2"+ 








K O* 

2^+ 


la 






















la + 2i Ja 























la - 2i Ja 











!"+ 











lb 











2"+ 











lb 























lb + 2i Jj 
























(34) 











lb - 2i Jfe 



(35) 



and 



where 



5. 



pm 



54 



Bo 
Bo 



(36) 



B^ 



























(37) 



and the Bo are 4x4 null matrices. 
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In this case we can find analytical solutions for the eigenvalues of Apm, which we 
write in terms of and /5+ as 

■^1,2 = o 



A 



3,4 



(38) 



^5,6 = 7a ± 2iJa, 

A7,8 = 76 ± 2iJ6, 
where T = 7^ + 7;,. 




Figure 5. The output spectral correlation S'°"*(Xp) + S°'"'*{Ym) of the low frequency 
modes, for k = 0.01, ja = 1, 7b = 2, Aa = Ja = 107a, and At = Jt = 27a. 

There are two types of instability which we may expect in this system and which 
would invalidate a linearised fluctuation analysis. The first is when one or more 
of the eigenvalues above has a negative real part and the second is the self-pulsing 
regime [351 EHl ET], which begins where there exist complex conjugate eigenvalues with 
real part equal to zero. Examining the expressions of (!38|) . we see that the last four 
can never cause any problems, but that some of the others could develop negative real 
parts. While analytical solutions for these in terms of the pump strength rather than 
the cavity field values can be found, these are extremely unwieldy. However, numerical 
analysis shows that the system is stable up to the same critical pumping as given in 
(p8|) . so we will give results in this regime. 



Bright entanglement 



15 




Figure 6. The output spectral correlation S°^'^{Xp) + S'°"*(F,„) of the high frequency 
modes, for k = 0.01, 7a = 1, 7;, = 2, = Ja = 107a, and Af, = Jb = 2^a- 



In terms of the quadrature variances used in section |H we now define 

Xm = Am + v4+ = Xi - X2, 

Y, = -I (Ap - A+) = Fi + Fa, 

Ym = -I [Am - A+) = Fi - F2, 

and similarly for the second harmonic quadratures, so that we can give expressions 
for the output spectral variances of these new quadratures. Because of the way the 
variables are defined, the easiest correlations to extract are those of ( fT2|) . which can be 
constructed from the variances of the four quadratures defined above, with 

V{X^ - X2) + V{Y^ + Fs) = V{Xm) + ^(Fp), 
y(Xi + X2) + - F2) = n^p) + V^(FJ. 

We find that, at least in the parameter regimes we have investigated numerically, it is 
the second of these correlations which violate the inequalities, which is different to the 
situation with, for example, coupled downconverters [2]. 

We show spectral results obtained numerically as the pumping strength is varied 
up to its critical value, and for different values of the cavity loss rates. In figure [5] and 
figure [6] we present results for the loss rate at the harmonic frequency being twice that 
of the fundamental, showing that the violations increase as the pumping increases up 
to the critical value. In figure [5] we see that the spectrum bifurcates and note that 
this could be advantageous experimentally as the region around zero frequency is often 
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Figure 7. The output spectral correlation S°"*{Xp) + S'°"*(Fm) of the low frequency 
modes, for k = 0.01, ja = 1, 7b = 7a/2, Aq = Ja = lO^a, and Af, = Jb = 2^a- 

swamped by technical noise. In figure [7] and figure M we give the spectra for the case 
where the loss rate at the fundamental frequency is twice that of the harmonic. The 
main difference is that the violations of the inequalities happen over a narrower range of 
frequencies, with both violations again increasing as the pump increases. In comparison 
with coupled downconversion below threshold |2], the entanglement produced here is 
with reasonably intense fields, as shown in figure [9], which could be a real advantage in 
some applications. 

7. Conclusions 

We have analysed the system of intracavity evanescently coupled second harmonic 
generation in terms of phase sensitive correlations which give evidence of continuous 
variable entanglement between different modes of the optical field. We have shown how 
matching the cavity detunings to the evanescent coupling rates fixes the quadrature 
angles for which the best violations of the inequalities occur. We see that below the 
self-pulsing threshold, the system exhibits a wide range of behaviour depending on the 
relative strengths of the pumps, the cavity loss rates, the detunings and the evanescent 
couplings. As all of these are experimentally tuneable, the device may be of use for 
applications which require varying degrees of entanglement to be available at different 
intensities, different frequencies and different quadrature phase angles. The entangled 
beams exit the cavity at different spatial locations and do not have to be separated before 
measurements can be made. As the system also produces single mode squeezing and can 
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Figure 8. The output spectral correlation S°^'^{Xp) + S'°"*(F,„) of the high frequency 
modes, for k = 0.01, 7a = 1,7;,= 7a/2, = Ja = 107a: aud Ai, = Jb = 2-/a- 




Figure 9. The intracavity intensities for each mode as a function of the ratio e/cc. 
The upper pair and the lower pair are for different ratios of the cavity loss rates at 
each frequency, with k = 0.01, 7a = 1, Aa = Ja = 107a, and Ab = Jb = 27a. 
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be built using integrated optics, it may prove to be more robust and useful than devices 
which rely on the relative stability and positioning of individual optical components. 
Finally, the fact that the entanglement is present with reasonably intense fields may 
prove to be a real advantage over devices based on nondegenerate downconversion, 
which experience phase diffusion in the region where the fields develop macroscopic 
intensities. 
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